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ABSTRACT 

Recently, a new planet candidate was discovered on direct images around the young (10-17 Myr) A-type star 
HD 95086. The strong infrared excess of the system indicates that, similar to HR8799, ft Pic, and Fomalhaut, 
the star harbors a circumstellar disk. Aiming to study the structure and gas content of the HD 95086 disk, and to 
investigate its possible interaction with the newly discovered planet, here we present new optical, infrared, and 
millimeter observations. We detected no CO emission, excluding the possibility of an evolved gaseous primordial 
disk. Simple blackbody modeling of the spectral energy distribution suggests the presence of two spatially separate 
dust belts at radial distances of 6 and 64 AU. Our resolved images obtained with the Herschel Space Observatory 
reveal a characteristic disk size of ~6'/0 x 5'.' 4 (540 x 490 AU) and disk inclination of ~25°. Assuming the same 
inclination for the planet candidate’s orbit, its reprojected radial distance from the star is 62 AU, very close to the 
blackbody radius of the outer cold dust ring. The structure of the planetary system at HD 95086 resembles the 
one around HR8799. Both systems harbor a warm inner dust belt and a broad colder outer disk and giant planet(s) 
between the two dusty regions. Modeling implies that the candidate planet can dynamically excite the motion of 
planetesimals even out to 270 AU via their secular perturbation if its orbital eccentricity is larger than about 0.4. 
Our analysis adds a new example to the three known systems where directly imaged planet(s) and debris disks 
coexist. 
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1. INTRODUCTION 

Recently, several new massive exoplanets at large orbital 
radii were directly imaged, e.g., around HR8799 (Marois et al. 
2008), /3 Pic (Lagrange et al. 2009), and Fomalhaut (Kalas et al. 
2008). They are the largest products of the planet formation 
process, which also produces many smaller planetesimals. 
While these smaller bodies cannot be detected directly, debris 
dust arising from their erosion can be observed via scattered 
light and/or thermal emission of the grains. All three examples 
mentioned harbor luminous debris disks. Giant planet(s) and 
the planetesimal belt/debris disk can interact with each other in 
several ways. A giant planet can sculpt the structure of the debris 
disk (Ertel et al. 2012), while the planetesimal disk can also 
influence the planet’s orbital evolution and long-term stability 
(Moore & Quillen 2013). 

Rameau et al. (2013) discovered a new planet candidate 
around HD 95086, at a projected separation of 56 AU. Its 
estimated mass of 4-5 Mj up makes it the lowest mass planet 
detected by direct imaging. The host star, similar to HR8799 
and P Pic, is an A-type star, belongs to a young association 
(Lower Centaurus Crux, de Zeeuw et al. 1999, hereafter LCC), 


13 ESA fellow. 


and exhibits prominent infrared excess indicative of a dusty 
circumstellar disk (Rhee et al. 2007; Chen et al. 2012; Rizzuto 
et al. 2012). The age estimates of the LCC association range 
between 10 Myr (Song et al. 2012) and 17 Myr (Pecaut et al. 
2012). While the disk of HD 95086 probably contains secondary 
dust, considering its young age, it is possible that the disk is an 
evolved gaseous primordial disk. Here, we study the structure 
and gas content of the HD 95086 disk and investigate its possible 
interaction with the newly discovered planet. Our analysis adds a 
new example to the three known systems where directly imaged 
planet(s) and debris disks coexist. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed HD 95086 with the Herschel Space Observa- 
tory (Pilbratt et al. 2010) using the Photodetector Array Camera 
and Spectrometer (PACS; Poglitsch et al. 2010), and the Spec- 
tral and Photometric Imaging Receiver (SPIRE; Griffin et al. 
2010). PACS observations were performed on 2011 August 1 
in mini- scan map mode with a medium scan speed (20" s _1 ). 
Measurements at scan angles of 70° and 110° were made both 
at 70 /z m and 100 /zm and this setup also provided four obser- 
vations at 160 /zm. Each measurement included four repetitive 
scan maps. Data processing was carried out with the Herschel 
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Figure 1 . Upper panels: Herschel PACS 70/100/160 /xm images of HD 95086. Lower panels: radial brightness profiles for HD 95086 and reference PSF profiles at 
the same wavelengths (see Section 3.2). 

(A color version of this figure is available in the online journal.) 


Interactive Processing Environment (HIPE; Ott 2010), v.9.2, 
using the pipeline optimized for bright sources. 1 // noise was 
removed by highpass filtering, after our target was masked to 
avoid flux loss. We used second-level deglitching to remove 
glitches. Mosaics (Figure 1) were created with pixel sizes of 
1". 1, 1M, and 2"1 at 70, 100, and 160 /zm, by combining the 
individual scan maps with a weighted average. SPIRE maps 
were obtained on 2011 August 16 at 250, 350, and 500 /z m in 
the small scan map mode with a repetition factor of two. Data 
reduction was performed with HIPE v9.2 using the standard 
pipeline script. The beam sizes were 1 8" 1 , 25 / . / 2, and 36'! 6 at 
250, 350, and 500 /z m, respectively, and the maps were resam- 
pled to pixel sizes of 6 ", 10", and 14". 

HD 95086 was observed with the InfraRed Spectro- 
graph (IRS) onboard the Spitzer Space Telescope on 2004 
February 4. Small 2x3 maps were taken using the low- 
resolution IRS modules, covering the 5.2-38 pm wavelength 
range with a spectral resolution of R = 60-120. We down- 
loaded the data processed with the pipeline version SI 8. 18.0 
from the archive and further processed them with the Spitzer 
IRS Custom Extraction Software (SPICE v2.5.0). We treated 
the data as a nodding measurement by using the two central 
map positions and subtracting them from each other. Then, we 
extracted the positive signal from a wavelength-dependent ta- 
pered aperture and averaged them. The final spectrum is plotted 
in Figure 2(a). 

On 201 1 April 22, we observed HD 95086 at the 345.796 GHz 
12 CO J = 3-2 line using the SHeFI/APEX2 receiver (Vassilev 
et al. 2008) mounted at the 12 m APEX telescope 14 (Giisten 
et al. 2006) (M-087.F-0001 program, PI: Th. Henning). We used 


14 This publication is based on data acquired with the Atacama Pathfinder 
Experiment (APEX). APEX is a collaboration between the 
Max-Planck-Institut fur Radioastronomie, the ESO, and the Onsala Space 
Observatory. 


the Fast Fourier Transform Spectrometer with 2048 channels, 
providing a velocity resolution of 0.42 km s _1 . An on-off ob- 
serving pattern was utilized with beam switching. The total on- 
source integration time was 29.4 minutes. The data reduction 
was performed using GILD AS /CL ASS. 15 For the final spec- 
trum, we subtracted a linear baseline from each individual scan 
and averaged them after omitting the noisy scans. 

We obtained a high-resolution optical spectrum of HD 95086 
with the Fiber-fed Extended Range Optical Spectrograph 
(FEROS; Kaufer et al. 1999) mounted at the 2.2 m MPG/ESO 
telescope in La Silla, Chile, on 2011 April 17. This instrument 
covers the wavelength range between 3500 and 9200 A in 39 
echelle orders with a resolution of R ~ 48,000. We used the 
“object-sky” mode, with one fiber positioned at the target and 
the other one on the sky. The integration time was 120 s. Data 
reduction, including bias subtraction, flat-field correction, back- 
ground subtraction, the definition and extraction of orders, and 
wavelength calibration, was performed using the FEROS data 
reduction system pipeline at the telescope. 

3. RESULTS AND ANALYSIS 
3.1. Stellar Properties 

To estimate the basic properties of HD 95086 and provide 
photospheric flux predictions at relevant mid- and far-IR wave- 
lengths, we modeled the stellar photosphere by fitting the opti- 
cal and infrared photometric data with an ATLAS9 atmosphere 
model (Castelli & Kurucz 2003). We used photometry from 
Hipparcos (Perryman et al. 1997), TYCH02 (H0g et al. 2000), 
and Two Micron All Sky Survey (Skrutskie et al. 2006) that 
were supplemented by W l band (centered on 3.4 /zm) photom- 
etry from WISE (Wright et al. 2010). Assuming that the inter- 
stellar reddening is negligible and adopting solar metallicity and 

15 http://iram.fr/IRAMFR/GILDAS/ 
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Figure 2. Left: color-corrected SED of the source and fitted models. Right: absolute V band magnitude vs. V — I c color diagram for known members of LCC with 
parallax error <10% (gray dots) and for HD 95086 (black square). Photometric data were taken from the Hipparcos catalog and dereddened using extinction values 
from Chen et al. (2011, 2012). 

(A color version of this figure is available in the online journal.) 


Table 1 

Disk Properties Derived from Quadratic Deconvolution (See Section 3.2) 



Original FWHM of the 
Source Before Deconvolution 

Disk Size After 
Quadratic Deconvolution 

Position 

Angle 

Incl. 


Major x Minor Axis (") 

Major x Minor Axis (") 

Major x Minor Axis (AU) 

(°) 

(°) 

PACS 70 jum 
PACS 100 /zm 
PACS 160 /zm 

8.1 ±0.1 x 7.7 ±0.1 

9.1 ±0.1 x 8.6 ±0.1 
12.5 ±0.5 x 12.2 ±0.4 

5.9 ±0.1 x 5.4 ±0.2 
6.2 ±0.2 x 5.4 ±0.2 
6.7 ± 1.0 x 6.1 ±0.9 

-530 x 490 
-560 x 490 
-600 x 550 

99 ±5 
1 13 ± 7 

23 ±5 
29 ±4 


log g = 4.0, our x 2 minimization yielded T Q ff = 7500zbl50K 
and a luminosity of Lt> 0 i = 7.0 zb 0.6L© (using a Hipparcos- 
based distance of 90.4 pc; van Leeuwen 2007). The derived T e ff 
is identical to that obtained by Chen et al. (2012) from color 
indices. 

To obtain log g , and v sin i from our FEROS spectrum, 
we used an iterative fitting method described in e.g., Szabo et al. 
(2011). First, we varied these three parameters using spectra 
from Munari et al. (2005) in the 4000-6200 A wavelength range, 
excluding the Hp region and Na D lines. Radial velocity was 
calculated from the cross-correlation function, convolving the 
measured spectrum with the best-fit Munari template, but setting 
log g to 1.0 and v sin i to 0. This yielded 17 zb 2 km s _1 . After 
transforming the measured spectrum to the laboratory system, 
and keeping T^f , log g , and v sin i in the global minimum, we 
fitted [Fe/H] in the second step. T e ff and log g was refitted in the 
third iteration. Our best-fit parameters are T ef f = 7750 dz 250 K, 
log g = 4.0 =b 0.5, v sin i = 20 zb 10 km s _1 , [Fe/H] = 
—0.25 zb 0.5. This T ci j supports our result obtained from the 
photometry (see above). Rameau et al. (2013) estimated a stellar 
mass of 1.6 M© for HD 95086. 

3.2. Herschel Maps 

HD 95086 appeared as an isolated bright source on our 
PACS images (Figure 1, upper panels). To estimate the disk 
size, we constructed point-spread functions (PSFs) using PACS 
observations of two calibrator stars ( a Boo, a Tau) without 
known infrared excess, measured and processed in the same 
way as HD 95086. The PSF maps were rotated to match the 


roll angle of the telescope at the time of observing HD 95086. 
The FWHM point- source sizes, as measured on these maps, are 
5" 5 zb an, 6"7 zb ai, and l(T/6 zb a 3 at 70, 100, and 160 /zm, 
respectively. We derived azimuthally averaged radial brightness 
profiles and averaged the individual measurements to obtain a 
reference PSF. By comparing the azimuthally averaged radial 
brightness profiles of HD 95086 with these reference PSFs, we 
found that the disk around HD 95086 is spatially extended at all 
PACS wavelengths (Figure 1, lower panels). 

We fitted elliptical Gaussians to the images of the source 
in order to derive the FWHM sizes of the major and minor 
axes and the position angles (see Table 1). Then, we estimated 
the disk size by quadratic deconvolution with the PSF, and 
calculated the inclination (at 70 and 100 /zm) assuming that the 
intrinsic shape/ structure of the disk is azimuthally symmetric. 
The contribution of stellar photosphere to the total flux was 
negligible (< 0.3%) in all bands, and thus we did not subtract the 
model photospheric emissions from the source images before 
the Gaussian fitting. Table 1 lists the final parameters and their 
uncertainties. The coarse spatial resolution at 160 /zm prevented 
us from determining a reliable position angle and inclination. 
Our results outline a large disk with a characteristic size of 
~6"0 x 5'.'4 (-540 x 490 AU). 

We measured the flux of HD 95086 on the individual 
PACS scan maps using an aperture radius of 18" and a sky 
annulus between 40" and 50". We calculated the average and 
RMS of the individual flux values, applied aperture correction, 
and calculated the final uncertainties of the photometry by 
quadratically adding the measurement errors and an absolute 
calibration uncertainty of 7% (Balog et al., submitted). Our 
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Table 2 

Measured and Predicted Fluxes 


Wavelength 

(/zm) 

Measured Flux Density 
(mJy) 

Instrument 

Predicted Photo spheric Flux Density 
(mJy) 

Reference 

3.35 

626.7 ± 24.7 

WISE 

592.3 

Wright et al. (2010) 

4.60 

336.8 ± 11.5 

WISE 

329.4 

Wright et al. (2010) 

9.03 

95.7 ± 13.0 

IRS 

90.7 

This work 

11.02 

67.3 ± 7.0 

IRS 

61.3 

This work 

11.56 

62.7 ± 3.0 

WISE 

55.8 

Wright et al. (2010) 

13.02 

51.2 ±5.4 

IRS 

44.2 

This work 

14.87 

44.7 ± 4.6 

IRS 

34.0 

This work 

16.99 

40.0 ± 4.4 

IRS 

26.1 

This work 

19.02 

43.1 ±5.0 

IRS 

20.9 

This work 

21.31 

46.7 ± 7.3 

IRS 

16.7 

This work 

22.09 

51.6 ±3.3 

WISE 

15.5 

Wright et al. (2010) 

23.67 

45.6 ± 2.0 

MIPS 

13.5 

Chen et al. (2012) 

24.48 

59.0 ± 6.0 

IRS 

12.6 

This work 

27.45 

77.3 ± 7.6 

IRS 

10.1 

This work 

30.50 

96.7 ± 10.0 

IRS 

8.1 

This work 

33.55 

137.9 ± 18.6 

IRS 

6.7 

This work 

60.00 

601.0 ±48.1 

IRAS 

2.1 

Moshir (1989) 

70.00 a 

690.1 ±48.6 

PACS 

1.5 

This work 

71.42 

655.0 ± 44.4 

MIPS 

1.5 

Chen et al. (2012) 

100.00 a 

675.1 ±47.6 

PACS 

0.7 

This work 

160.00 a 

462.4 ± 32.7 

PACS 

0.28 

This work 

250.00 

213.4 ± 12.9 

SPIRE 

0.12 

This work 

350.00 

120.3 ± 8.7 

SPIRE 

0.06 

This work 

500.00 

63.6 ± 10.2 

SPIRE 

0.03 

This work 


Note. a The emission is spatially resolved at these wavelengths. 


target was clearly detected as a point source at all SPIRE 
wavelengths. We performed photometry with an aperture radius 
of 22", 30", and 42" at 250, 350, and 500 \i m. Background 
levels were estimated in annuli extending from 60" to 90". 
The final uncertainties were derived as the quadratic sum of 
the measurement errors and the overall calibration uncertainty 
of 5.5% for the SPIRE photometer (Bendo et al. 2013). Our 
Herschel photometry for HD 95086 is listed in Table 2. 

3.3. Modeling of Dust Distribution 

We compiled the spectral energy distribution (SED) of 
HD 95086 by combining the new PACS and SPIRE fluxes 
with infrared photometry from the literature. For the fitting 
process the IRS spectrum was sampled in 11 bins. Table 2 
summarizes the collected IR data. The resulting SED with 
the stellar photosphere model (Section 3.1) is presented in 
Figure 2(a). The IRS spectrum shows that IR excess is present 
even at ~10/rm. No silicate features are apparent, implying 
the depletion of small grains. We assumed optically thin dust 
emission in the modeling, and, as a first simple model, we fitted 
the excess by a single temperature modified blackbody, where 
the emissivity is equal to 1 at k < Xq and varies as (A/Ao) - ^ 
at A. > A.Q- Following Williams & Andrews (2006), we adopted 
A 0 = 100 fi m. We used a Levenberg-Marquardt algorithm in 
the modeling, and utilized an iterative method to compute color 
corrections (Moor et al. 2006). Table 3 lists the best-fit pa- 
rameters and the reduced x 2 of the fitting. Figure 2(a) shows 
that the best-fit model underestimates the observed excess, espe- 
cially at shorter wavelengths. This discrepancy can be explained 
if the emitting dust grains — similar to many other debris sys- 
tems (e.g., Morales et al. 2011) and our solar system — are dis- 
tributed in two spatially separated rings. With this assumption, 
we used a two-component model, where grains in the warmer 


Table 3 

Disk Properties Derived from SED Fitting 



Single Temperature 

Two Temperature 



Warm Comp. 

Cold Comp. 

Xr 

4.6 

0.9 


Tdust (K) 

64 ± 1 

187 ± 26 

57 ± 1.5 

p 

0.17 ±0.06 

0 

0.41 ± 0.07 

r dust (AU) 

50 ± 2.4 

5.9 ± 1.6 

63.7 ± 4.4 

/dust (10 3 ) 

1.5 ±0.1 

0.14 ±0.10 

1.4 ±0.2 


component act like blackbodies, while the emission of the outer 
ring is described by the modified blackbody as defined above. 
As Figure 2(a) and the improved x 2 value (Table 3) indicate, 
this model fits the SED better over the whole studied wavelength 
range. The fractional luminosity of the rings was computed as 
/dust = ^dust/^boi- The radius of the dust ring(s) was estimated 
using the following formula (Backman & Paresce 1993): 

^dust = / LboA 0 ' 5 / 278K\ 2 

AU ( L q ) ( Tdust ) ' ' 


Because of the blackbody assumption, the resulting rd us t values 
are lower limits. The obtained fundamental disk properties are 
listed in Table 3. 

The dust mass of the disk was estimated using the following 
formula: 


3Tdust — 


F v d 2 

Bv(7dust>v 


( 2 ) 


where F v is the measured SPIRE flux at 500 /im, d is the 
distance to the source (90.4 pc), k v = /co(v/vo)^ is the mass 
absorption coefficient, and B v = ((2v 2 £7d us t)/c 2 ) is the Planck 
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function using the Rayleigh-Jeans approximation. Assuming a 
/c 0 = 2 cm 2 g -1 at v 0 = 345 GHz (e.g., Nilsson et al. 2010) and 
taking j3 = 0.41 and r dust = 57 K from Table 3, we derived a 
dust mass of 0.5 ±0.1 M®. 

3.4. CO Data 

No CO emission was detected in our APEX spectrum. We 
converted antenna temperatures to flux densities using conver- 
sion factors from the APEX Web site 16 and computed an upper 
limit for the CO(3-2) line flux (*Sco(3-2)) as AmsAi ’ a/A) where 
/ rms = 0.95 Jy is the measured rms noise at the systemic veloc- 
ity of the star, Av = 0.42 km s _1 is the velocity channel width, 
and N is the number of velocity channels over an interval of 
10 km s _1 , which covers the expected line width. The resulting 
upper limit is 1.95 Jy km s _1 for Sco( 3 - 2 )- Assuming optically 
thin emission and using a range of excitation temperatures be- 
tween 20 K and 60 K (the latter value corresponds to the dust 
temperature in the outer ring; see Section 3.3), the upper limit 
for the total CO mass is between 1 .4 x 10 -4 and 1 .7 x 10 -4 M® . 

4. DISCUSSION 

Analyzing Hipparcos astrometric data, de Zeeuw et al. 
(1999) identified a total of 180 probable members in the LCC 
association. Among their proposed members, HD 95086 has 
the lowest membership probability of 41%. In order to re- 
check the membership status of HD 95086, we used our new 
radial velocity and Hipparcos astrometric data to compute a 
Galactic space motion of U = —10.5 ± 1.0, V = —22.3 ± 1.8, 
W = —5.2 ± 0.4 km s _1 . According to Chen et al. (2011), 
the characteristic space motion of LCC is U = —7.8 ±0.5, 
V = — 20. 7 ±0.6, W = — 6.0 ±0.3 km s _1 , thus there is a good 
kinematic match between our target and LCC. The position of 
HD 95086 in the color-magnitude diagram in Figure 2(b) fits 
well to the locus defined by LCC members, indicating that its 
isochrone age is consistent with that of the cluster. Based on 
these findings, we propose that HD 95086 is a very probable 
member of LCC. With its distance of 90.4 pc, HD 95086 is one 
of the closest members of the association. 

Pecaut et al. (2012) derived an age of 17 Myr for LCC based 
on isochrone fitting while Song et al. (2012) proposed 10 Myr 
from lithium content of late-type LCC stars. The latter estimate 
is close to the timescale on which primordial gas is believed 
to disappear (e.g., Mamajek 2009). Thus, if the 10 Myr age 
estimate proves to be true, we cannot exclude that this disk 
would contain substantial amount of primordial gas similar to 
transitional disks. Converting our upper limit on the CO mass to 
total gas mass assuming a canonical H 2 /CO abundance ratio of 
10 4 , the estimated gas mass in the disk of HD 95086 is <0. 12 M® 
(<3.7 x 10 -6 M q ). This is several orders of magnitude below the 
typical Herbig Ae and transitional disk masses. This result and 
the low gas-to-dust mass ratio of <0.25 implies that HD 95086 
harbors a gas-poor debris disk. 

The debris disk of HD 95086 has a high fractional luminosity 
of 1.5 x 10 -3 . With this value, it belongs to the top 10 highest 
fractional luminosity debris systems in the solar neighbourhood 
(< d < 120 pc; Moor et al. 2006; Rhee et al. 2007). Interestingly, 
these systems have other common properties: they are younger 
than 100 Myr and most of them belong to different young mov- 
ing groups and associations. We found that the HD 95086 disk 


16 http://www.apex-telescope.org/telescope/efficiency/ 


contains two dust belts: a warm one with a characteristic tem- 
perature of 187 K, and a colder one with 57 K. Such structures 
seem to be relatively common. Morales et al. (201 1) found that 
a significant fraction of debris disks harbors two spatially dis- 
tinct dust components. The characteristic dust temperatures for 
the inner and outer dust belts (190 K and 60 K) were found 
to be very similar independently of the host stars’ luminos- 
ity, implying that the formation of dust belts is controlled by a 
temperature-sensitive mechanism. Several other members of the 
young massive debris disk group also share similar properties 
(e.g., Su et al. 2009; Roberge et al. 2013; Donaldson et al. 2013). 
Because the characteristic temperature of the warm component 
is slightly above the ice evaporation temperature, Morales et al. 
(2011) proposed that the formation of the warm component is 
related to sublimation of icy planetesimals crossing the snow 
line, or due to collisions in an asteroid belt-like system formed 
just interior to the snow line (see also Martin & Livio 2013). One 
of these mechanisms could also work in the case of HD 95086. 

To determine the relative location of the two dust belts and 
the planet candidate, we need to constrain the inclination of the 
planetary orbit. From our resolved Herschel observations, the 
disk inclination is ~25°. The star’s equatorial plane is probably 
not very different from this, because the observed v sin i = 
20 km s _1 is unusually low compared to stars with similar 
masses (Zorec & Royer 2012), implying a low stellar inclination. 
Assuming that the disk and the planet’s orbit are coplanar, and 
taking the projected orbital radius of 56 AU from Rameau et al. 
(2013), the deprojected orbital radius would be 62 AU. We 
derived a radius of 64 AU for the outer dust belt. Since we 
adopted blackbody grains, this is a lower limit, because smaller 
grains with the same temperature can be located farther from the 
star. Thus, it is possible that the planet candidate is situated just 
inside the cold outer dust belt. Further monitoring of the planet 
candidate will reveal whether it orbits within the outer dust belt. 
If it stays inside, it might sculpt the inner edge of the belt, and 
may induce azimuthal asymmetries in the dust distribution. We 
note that the planetary system of HD 95086 resembles the one 
around HR8799. Both systems harbor a warm inner dust belt 
and a broad colder outer disk and giant planet(s) between the 
two dusty regions. 

The short-lived dust grains in debris disks are believed to be 
replenished by collisions between larger bodies. For destructive 
collisions, the collision velocity must exceed a critical value and 
it requires a dynamical perturbing force. According to theory, 
this perturbation is linked to 1000 km-sized planetesimals 
formed within the planetesimal ring (self- stirring; Kenyon & 
Bromley 2008), or the presence of planet(s) in the system 
(planetary stirring; Mustill & Wyatt 2009). The formation of 
1000 km sized planetesimals is a slow process at large orbital 
radii. Adopting an age of 17 Myr for the system and using 
Equations (27) and (41) from Kenyon & Bromley (2008), we 
found that, even in a disk with an unusually high initial surface 
density (10 times higher than the minimum-mass solar nebula), 
the formation of such large bodies is limited to <70 AU. 
Our PACS images show that the disk around HD 95086 is 
spatially extended with a diameter of 540 AU. While large 
dust grains have the same spatial distribution as the parent 
planetesimals, the distribution of small grains (smaller than the 
blowout limit of ~ 1.8 /rm, estimated following Wyatt 2008) can 
be more extended, because they are expelled by stellar radiation 
pressure and form a dust halo. However, the contribution of such 
small grains to the 100 and 160/xm flux is small, suggesting 
that HD 95086 harbors an extended planetesimal belt whose 
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dynamical stirring cannot be explained by self- stirring. Mustill 
& Wyatt (2009) claim that planetary stirring can be faster 
than self-stirring in the outer regions of planetary systems. 
Using their formulae, we found that the candidate planet can 
dynamically excite the motion of planetesimals even at 270 AU 
via their secular perturbation if its orbital eccentricity is larger 
than ~0.4. Alternatively, without effective stirring, collisions 
with low velocity can lead to merging grains, thus we cannot 
exclude that the outer disk region contains primordial dust 
as well. 
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